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Abstract 

We use Omnes representations of the form factors /+ and /o for exclusive semileptonic 
B — > % decays, paying special attention to the treatment of the B* pole and its effect 
on /+. We apply them to combine experimental partial branching fraction information 
with theoretical calculations of both form factors to extract \V u b\. The precision we 
achieve is competitive with the inclusive determination and we do not find a significant 
discrepancy between our result, \V u b\ = (3.90±0.32±0.18) x 10 -3 , and the inclusive 
world average value, (4.45 ±0.20 ±0.26) x 10~ 3 [1]. 

1 Introduction 

The magnitude of the Cabibbo-Kobayashi-Maskawa matrix element V u b can be determined from 
both inclusive and exclusive semileptonic B meson decays. There has been a recent dramatic 
improvement in the quality of the experimental data for the exclusive decays [2-6], coupled with 
the appearance of the first dynamical lattice QCD and improved lightcone sumrule calculations 
of the relevant form factors [7-12]. Dispersive approaches were combined with lattice results 
in [13] and with leading order heavy meson chiral perturbation theory and perturbative QCD inputs 
in [14]. The appearance of the first partial branching fraction measurements for B — > nlv [2] 
made it possible [15] to combine dispersive constraints with experimental differential decay rate 
information and theoretical calculations of both form factors in limited regions of q 2 in order to 
improve the determination of \V u b\. In [16] it was shown that the quality of the inputs now makes 
it possible for the exclusive determination to compete in precision with the inclusive oneQ. Thus 
the compatibility of the two determinations becomes an interesting issue. 

To perform the exclusive \V u b\ extraction one needs a model-independent parametrisation of the 
form factors. In [16] a parametrisation inspired by dispersive bounds calculations was used. An 
alternative simple parametrisation using a multiply-subtracted Omnes representation for / + , based 
on unitarity and analyticity properties, was employed in [18]. A shortcoming in the treatment of the 
B* was pointed out in [19]. In this letter we have addressed this by improving the treatment of the 
B* within the Omnes framework. We have also incorporated the scalar form factor /o in a simulta- 
neous analysis and examined the possible effects of correlations among lattice inputs. Finally, we 
have taken advantage of new experimental data from the BaBar 12-bin untagged analysis [6]. 

The outcome is that the precision achieved for \V u b\ is indeed competitive with the inclusive de- 
termination and that we do not find a significant discrepancy between our result, \V u b\ = (3.90 ± 
0.32 ±0.1 8) x 10~ 3 , and the inclusive world average value, (4.45 ±0.20 ±0.26) x 10~ 3 [1]. 

^ee [17] for updates of the fits in [16]. 
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2 Omnes Parametrisations 



In our previous work [18, 20] with the Omnes parametrisation [21, 22] for the form factor f+(q 2 ), 
we treated the B* as a bound state and took the B% elastic scattering phase shift to be % at threshold, 
Sth = ( m B + tnx) 2 - By using multiple subtractions and approximating the phase shift by n from 
to infinity, this led to a parametrisation: 

f + (q 2 ) = — i-jff [M Si )( S& - Si )] ai ^\ (1) 

with n+l subtractions at q 2 G {^o, *i , • • • ,s n }, below threshold (the CCi(q 2 ) are defined in equa- 
tion © below). This parametrisation requires as input only the form factor values {f+(qf)} at 
n+l positions qf. 

Using this parametrisation in a combined fit to experimental data and theoretical form-factor cal- 
culations (lattice QCD and lightcone sumrules) allows an extraction of \V u b\ with precision com- 
petitive to the inclusive determination. This parametrisation and others were compared in refer- 
ence [19] where the form factor /+ was determined by fitting BaBar experimental partial branching 
fraction data in 12 bins [5,23] and using \V u b\ determined from Unitarity Triangle fits. Good agree- 
ment was found between the Omnes parametrisation of equation (OQ) and parametrisations using 

for two choices of to. The coefficients a n satisfy the dispersive constraint Y,n a n — 1 H^]. Ex- 
pressions for P and can be found in [16]. When we use the parametrisation in equation ©, we 



will set to = Sth(l — v 1 — ^maxAth)* which is the 'preferred choice', labelled BGLa, in [19] (this 



choice for to ensures that \z\ < 0.3 for < q 2 < g^ax)- We will refer to the parametrisation using 
this functional form as the z-expansion or ZE below. 

Fits for f+(q 2 ) using the Omnes and ZE parametrisations deviated from each other by a few per 
cent only in the largest q 2 region, close to g^ax' which has little influence on the decay width and 
\Vub\ (see figure 2 in [19]). This is also the region where there is no theoretical information on 
the form factors. In figured] we show a similar comparison, including bands showing statistical 
fluctuations arising from the fits. We have fitted the same dataset as in [18], but replacing the 
5-bin BaBar untagged analysis [4] with the updated 12-bin results from [6]. The ZE fit has been 
performed truncating the power series in equation at n = 2, for comparison with figure 2 in [19], 
or n = 3, so that all fits have the same number of parameters. The green dashed lines show the 
Omnes fit using equation (0Q). The plot shows that once fluctuations are taken into account the 
differences are not significant. 

Nevertheless, we show here that by treating the B* explicitly as a pole of the form factor, we can 
understand and reduce the small deviation in the central fits at large q 2 . This is illustrated by the 
solid blue lines in figure [U We achieve this without altering the main results obtained for \V u b\ 
and / + in the q 2 region where theoretically calculated values lie. The new parametrisation, shown 
below in equation ©, is obtained from equation (OQ) by replacing s t h with m|*. As before, the 
parametrisation relies only on very general properties of analyticity and unitarity and so, although 
simple, is well-founded. 
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Figure 1 Comparison of fits to f+(q 2 ) using Omnes or ZE parametrisations. Each fit is plotted with its own 
error bands, but normalised by the central fit for the ZE. Thus, the horizontal line at 1 and the grey band 
show the ZE fit with its 68% statistical error band, the solid blue lines indicate the Omnes fit of equation ^} 
and the green dashed lines show the Omnes fit of equation £[]). The left-hand plot uses a ZE fit with three 
parameters and thus can be compared to figure 2 in [19], while the right-hand plot uses a four-parameter ZE 
fit. All Omnes fits use four parameters (four subtraction points). 



To obtain the new parametrisation we observe that, if f+(q 2 ) has a pole at q 2 = nig*, then J?(q 2 ) = 
(m|* — q 2 )f + (q 2 ) has no poles and satisfies 

-— — — = exp(2i5 1/2jl (s)), s>s th (3) 

where 8u is the phase-shift for elastic %B — > %B scattering in the isospin / and total angular mo- 
mentum J channel. This is because /+ satisfies a similar equation as required by Watson's theo- 
rem [24] and we have multiplied it by a real function. An (n+l)-subtracted Omnes representation 
can now be written for J^(g 2 ), with q 2 < s±, which reads: 



^{q 2 ) = [U[^(^r q } J exp {l s (q 2 ; s , . . .,s n ) Y\(q 2 - Sj )j, (4) 
I s (q ; SQ, ...,s n ) = - / r r (5) 

KJs th {s-So)---{S-S n ) S-q 1 

n s — s ■ n 
cci(s)= n — — -, cd(sj) = Sij, £a,-(s) = l. (6) 

j=o,#j s < s j i=o 

This representation requires as input the phase shift 5 1 /2,i(*) plus the values {J^"(s ( )} at n + 1 
positions {si} below the tzB threshold. For sufficiently many subtractions, we can approximate 
8\/2,\{s) by zero above threshold (see appendix lA)) . In this case we obtain, 

Mi 2 ) = -T^—zfl [M*i) ai{q2) • a) 

This amounts to finding an interpolating polynomial for \nJP(q 2 ) = ln[(m|* — q 2 )f+(q 2 )] pass- 
ing through the points &(si) = (m|* — Similarly, our earlier parametrisation in equa- 
tion (0Q) used an interpolating polynomial for ln[( s± — q 2 )f+ (q 2 )] ■ While one could always propose 
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a parametrisation using an interpolating polynomial for \n[g(q 2 )f + (q 2 )} for a suitable function 
g(q 2 ), the derivation using the Omnes representation shows that taking g(q 2 ) = m|* — q 2 (and 
equally s± — q 2 ) is physically motivated. From here onwards we will use equation CD) as our 
preferred parametrisation for f + . 

When using our parametrisation in the extraction of \V u b\, we make 4 subtractions. This is suf- 
ficient to justify using no information about the phase shift beyond its value at s±. To check 
this, we have put in a model for the Bit phase shift [20] and confirmed that induced changes 
in our results are much smaller than the fluctuations produced by the errors in our inputs. This 
can be understood because with four evenly-spaced subtractions at {0, 1/3, 2/3, l}<7 max > me factor 
exp [l§ x Yl'j^oil 2 ~ s j)] m equation © given by this model deviates from unity by no more than 
6 x 10 -4 for < q 2 < g max (and, of course, is unity at each subtraction point). 

Since /+ and and the scalar form factor fo satisfy the constraint /+(0) = /o(0) we will combine 
theoretical inputs for / + and fo with experimental B — > nlv partial branching fraction information 
to check the effect on the extracted value of \V u t\. We will investigate the effect of using the / + 
information alone or using both form factors. 

For the scalar form factor fo there is no pole below threshold, so that we will use an Omnes formula 
like equation © for fo(q 2 ), with & — > fo and Si^i ~~ > ^1/2,0- F° r sufficiently many subtractions, 
we can approximate Si no by zero above threshold. Our recent analysis of the scalar form fac- 
tor [25] for B — > % decays suggested the existence of a resonance with mass around 5.6 GeV. This 
could be incorporated in an Omnes parametrisation like that of equation ©, but (as we have con- 
firmed) has negligible effect on \V u b\ and /+ in our fit, producing only a small increase of around 
7% in the value of fo close to g max . 

3 Application to \V u b\ 

We have used experimental data for the partial branching fractions of B — > Ttlv decays in q 2 bins 
from both tagged and untagged analyses. The tagged analyses from CLEO [2], Belle [3] and 
BaBar [5] use three bins, while BaBar's untagged analysis [6] uses twelve. CLEO and BaBar 
combine results for neutral and charged 5-meson decays using isospin symmetry, while Belle 
quote separate values for 5° — > n~l + V[ and B + — > 7T°/ + V/. For our analysis, for the three-bin data, 
we have combined the Belle charged and neutral 5-meson results and subsequently combined 
these with the CLEO and BaBar results. The resulting input values can be found in table II of [18]. 
Since the systematic errors of the three-bin data are small compared to the statistical ones, we have 
ignored correlations in the systematic errors and combined errors in quadrature. For the 12-bin 
BaBar data [6], complete correlation matrices are available in the EPAPS database [26] for both 
statistical and systematic errors and we have used these in our fits (we used the results corrected 
for final state radiation effects). We have assumed no correlation between the untagged and the 
tagged analyses. 

When computing partial branching fractions, we have used T g o = 1/rxot = (1.527 ± 0.008) x 
1 12 s [ 1 ] for the B° lifetime. 

Since the effects of finite electron and muon masses are beyond current measurement precision, the 
experimental results provide information on the q 2 shape of / + . Theoretical calculations provide 
information on /+ and fo. 

We use the lightcone sumrule (LCSR) result /+(0) = / (0) = 0.258 ±0.031 [12] and lattice QCD 
results from dynamical simulations at larger q 2 from HPQCD [7] and FNAL-MILC [8-11]. The 
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FNAL-MILC results [8-11] are still preliminary. Therefore we use the three f+(q 2 ) values quoted 
in [16] and read off three values for fo(q 2 ) at the same q 2 points from [9]. These are 

/ (15. 87 GeV 2 ) =0.425 ±0.033 

/ ( 18.58 GeV 2 ) =0.506 ±0.037 (8) 

/ (24.09GeV 2 ) =0.800±0.067 

The errors shown are statistical. A further 11% systematic error should be added. 

We implement the fitting procedure described in [18] using four evenly-spaced Omnes subtraction 
points at {0, 1/3,2/3, 1}^^ (with ^-squared function given in equation (10) of [18]), with the 
obvious changes to incorporate /q. As before, we have assumed that the lattice input form factor 
data have independent statistical uncertainties and fully-correlated systematic errors. We have not 
assumed correlations between results for / + and /o, though we will comment further on this below. 
Furthermore, we ignore possible correlations between the HPQCD and FNAL-MILC lattice inputs. 
These correlations are unknown and we showed in [18] that unless they are very strong they will 
have little effect on \ V U £,\. 

The best-fit parameters are 

\V ub \ = (3.90 ± 0.32) x 10- 3 

/+(0)=/ (0) = 0.226±0.022 

/+(<7max/3) = 0.417±0.039 

f+(2q 2 m j3) = 0.941 ±0.064 

/+(<Zmax) = 7.29±1.28 w 

/o(?Lx/3) = 0.342 ±0.053 

fo(2q 2 ma j3) = 0.508 ±0.040 

/oteLJ = 1.09±0.21 

The fit has ^ 2 /dof = 0.62 for 28 degrees of freedom, while the Gaussian correlation matrix can be 
found in appendix [Bj 

In figure[2]we show the fitted form factors, the differential decay rate calculated from our fit and the 
quantities log[(m|* — q 2 )f+(q 2 )/mg t ] and P(j)f+ where the details of the fit and inputs can better 
be seen. The dashed magenta curve in the P(j)f+ plot is a cubic polynomial fit in z to the output 
from our analysis. We note that the sum of squares of the coefficients in this polynomial safely 
satisfies the dispersive constraint £„a 2 < 1 [16]. 

Compared to our previous results [18] we find that the central value of \V u b\ decreases by 3% 
compared with an error of around 8%. Similarly, the central values of /+(0) and f+iq^x) move 
up by around half their errors, while f+iq^x/i) increases by an amount comparable with its error. 
At 2g 2 nax /3, in the neighbourhood of which most of the form factor data is concentrated, there is 
hardly any change. The result for /o(<?max) agrees with that obtained in our recent analysis of the 
scalar form factor alone [25]. We make some remarks on these results: 

• We have checked that the changes in the results for /+(0), f+iq^ax/l) and \V u b\ stem from 
using the updated BaBar untagged data. 

• We have checked that the change in f+iq^x), which has little effect on the shape of the form 
factor in the q 2 range where experimental and theoretical information exists, arises from our 
use of the new Omnes parametrisation of equation CD) and reflects the existence of a pole in 
/+ Mq 2 = ml*. 
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Figure 2 Results obtained from the fit to experimental partial branching fraction data and theoretical form 
factor calculations. The top left plot shows the two form factors with their error bands, the lattice and LCSR 
input points (dots: green LCSR, red HPQCD, blue FNAL-MILC) and 'experimental' points (black triangles, 
upward-pointing for tagged and downward pointing for untagged data) constructed by plotting at the centre 
of each bin the constant form factor that would reproduce the partial branching fraction in that bin. The top 
right plot shows the differential decay rate together with the experimental inputs. The bottom plots provide 
more details of the inputs and fits by showing on the left log[(ra|» — q 2 )f+{q 2 )/m\ t ] as a function of q 2 , 
and on the right P<pf + as a function of — z. The dashed magenta curve in the bottom right plot is a cubic 
polynomial fit in z to the Omnes curve. 

• Since we do not know the correlations between the lattice input data we have also performed 
a fit neglecting all correlations in these inputs. We find that \V u b\ increases by an amount 
0.18 x 10 3 , which we will quote as a systematic error in our determination. We observe 
that knowledge of the correlations will be needed for more precise determinations of \V u t,\. 

• The inclusion of fo in the analysis has no visible effect in our results for / + and \V u b\. 
This is not surprising given that the number of input data affecting / + is much bigger than 
that affecting fo and that the parametrisation allows the data to determine each form factor 
independently apart from the constraint at q 2 = 0. The covariance matrix given in appendix|B] 
shows this freedom, having negligible correlations between / + and fo at q 2 ^ 0. Correlations 
linking /+ and fo in the lattice QCD inputs could modify the central values in © by an 
amount comparable to their errors as we have confirmed by fully-correlating the systematic 
errors between them. As an example, for \V u b\ we find a central value of 4.15 x 10~ 3 . Since 
we do not know the actual correlation informatiorH for the lattice data, we do not present 

2 It is reasonable to expect correlations not only in the systematic error but also in the statistical ones, since /+ and 
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these numbers. 

• Because of the freedom allowed by the Omnes parametrisation of /+ and fo, one may wonder 
whether or not heavy quark symmetry (HQS) relations between the form factors at q^ ax 
are satisfied. Some earlier parametrisations were explicitly constructed to satisfy the HQS 
scaling relation /+(<7m ax ) //o(?max) ~ m #> f° r example dipole/pole forms [27-29], and these 
have been widely used. From our fit we calculate 

= 1.3±0.4GeV _1 (10) 

fin- 
to be compared to the corresponding quantity in D — > % exclusive semileptonic decays, 1.4 ± 
O.lGeV 1 extracted from the unquenched lattice QCD results in [30]. This agreement is 
reassuring but our determination of the ratio in B — ► % decays has a further uncertainty of 
around 10% arising from our incomplete knowledge of the correlations in the lattice inputs. 

• Heavy quark effective theory in the soft-pion limit predicts [31], 

fo(m 2 B ) = f B /f n + 0{\/m 2 ) » 1.4(2) (11) 

where we have used f B = 189(27) MeV [32]. Our fit for /o(<?max) m equation © is compat- 
ible within errors. 

• Applying soft collinear effective theory (SCET) to B — > kk decays allows a factorisation 
result to be derived which leads to a model-independent extraction of the form factor (mul- 
tiplied by \V u b\) &t q 2 = [33]. We quote the result from our fit, 

|Kfc|/ + (0) = (8.8±0.8)xl0- 4 , (12) 

which compares well with \V u b\f + (0) = (7.6 ± 1.9) x 10 -4 quoted in [17]. This also agrees 
with the value |V^|/ + (0) = (9.1 ±0.7) x 10~ 4 [19] obtained by fixing \V ub \ from global 
CKM unitarity triangle fits and fitting to the BaBar 12-bin data [23]. 

• We noted above possible effects of correlations in the lattice data. Other sources of system- 
atic variation in the result for \V u j,\ arising from uncertainties in the theoretical form factor 
inputs at or near q 2 = were considered in [18] and shown to be safely covered by the 
statistical uncertainty. 

4 Conclusion 

We have updated our previous analysis of exclusive B — > % semileptonic decays, based on Omnes 
dispersion relations. The principal change is to improve the treatment of the B* and its effect on 
the form factor f + . We have also incorporated the scalar form factor fo in a simultaneous analysis 
and examined the possible effects of correlations among lattice inputs. Finally, we have taken 
advantage of new experimental data from the BaBar 12-bin untagged analysis [6]. We extract a 
value 

\V ub \ = (3.90±0.32±0.18) x 10~ 3 . (13) 

fy are linear combinations of temporal and spatial components of vector current matrix elements. 



J_ A^max) 
m B /o(^max) 
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The first error above is statistical arising from the chi-squared fit. The second is a systematic error 
to account for current partial knowledge of correlations in the lattice input data. The precision for 
| V u j,\ is comparable with that of the inclusive determination and we do not find a significant discrep- 
ancy between our result and the inclusive world average value, (4.45 ±0.20 ±0.26) x 10~ 3 [1]. 

Finally we would like to stress that the Omnes parametrisation is physically motivated and simple 
and provides a robust framework for a precise exclusive determination of \V u i,\. 
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A Choice of 8ij(s t ^) in the Omnes Representation 

In this appendix we provide more details on some aspects of the Omnes representation of the form 
factors. This builds on the discussion in the appendix of [34]. 

The scattering matrix T depends on exp(2?'<5) and thus one has the freedom to add factors of kn to 
the phase shift, for integer k, without modifying the T matrix. However, the Omnes representation 
of the form factor certainly depends on the specific value of k. Indeed adding kn to 8 leads to 

f n \ f n ^ ( W n (s s) a -i^ 2 ^\^ 

exp|Wxn(^ 2 -^)}=exp|/ 5 xn(^ 2 -^)j I J ^ 

which induces an unphysical kth order pole in the form factor at s±. 

Now consider <^(q 2 ) = (nig* — q 2 )f+(q 2 ), which has no poles in < q 2 < s t h- Its Omnes rep- 
resentation should not induce a pole at s± and therefore we should set k = in equation (fT4l) 
above. This is equivalent to setting 5i/2,i (*th) = 0. With enough subtractions, we can then take 
^1/2,1 ( s ) = inside the integral because only the region close to will be important, leading to 
the result presented in equation ©. 

This choice for 5j (^th) does not contradict Levinson's theorem, which fixes only the difference 

5i/2,i(°°) -5i/ 2 ,i(sth) = K{n z -n p ) (15) 

where n z (n p ) is the number of zeros (poles) of the scattering matrix T on the physical sheet. The 
usual convention [35] is to set 5^/2,1 (^th) = ^ n p an d ^1/2,1 (°°) = %n z- However, we use a different 
convention which follows from the discussion above on the effect of adding multiples of n to the 
phase shift. Our choice is b\ n\ (*th) = which therefore also implies that 5^2,1 (°°) = 0. 

In our previous work [18, 20, 34], we assumed that / + had no poles. With the usual convention 
for Levinson's theorem that 5i/2,i(sth) = ft, we developed a pole for / + at s t h which was not 
discarded since it mimicked the B* pole's effects on the form factor because m|* ta s t h- We already 
commented on this in the appendix of [34]. 
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B Correlation Matrix 



Here we give the correlation matrix of fitted parameters corresponding to the best-fit parameters in 
equation ©. 

/ 1 -0.43 -0.91 -0.81 -0.58 -0.04 0.00 
1 0.20 0.50 -0.04 0.10 0.00 
1 0.76 0.61 0.02 0.00 
1 0.36 0.05 0.00 
1 0.00 0.00 
1 0.32 
1 

V 
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